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Abstract 

Different approaches in quantifying environmentally-induced decoherence are considered. We 
identify a measure of decoherence, derived from the density matrix of the system of interest, 
that quantifies the environmentally induced error, i.e., deviation from the ideal isolated-system 
dynamics. This measure can be shown to have several useful features. Its behavior as a function of 
time has no dependence on the initial conditions, and is expected to be insensitive to the internal 
dynamical time scales of the system, thus only probing the decoherence-related time dependence. 
For a spin-boson model — a prototype of a qubit interacting with environment — we also demonstrate 
the property of additivity: in the regime of the onset of decoherence, the sum of the individual 
qubit error measures provides an estimate of the error for a several-qubit system, even if the qubits 
are entangled, as expected in quantum-computing applications. This makes it possible to estimate 
decoherence for several-qubits quantum computer gate designs for which explicit calculations are 
exceedingly difficult. 

PACS numbers: 03.67.Lx, 03.65.Yz 
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I. INTRODUCTION 

Dynamics of open quantum systems has increasingly attracted the at- 

tention of the community of scientists in diverse fields, working on realizations of quan- 
tum information processing. Recent interest in quantum computing has stimulated studies 



10 



11 



12 



13 



15 



16 



17 



19| of environmental effects that cause small 



deviations from the isolated-system quantum dynamics. To perform large-scale quantum 
computation, environment- induced relaxation/decoherence effects during each short time 
interval of "quantum-gate" functions must be kept below a certain threshold in order to 



allow fault-tolerant quantum error correction 
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26|,|27|. The reduced 



density matrix of the quantum system, with the environment traced out, is usually evalu- 

3, El- 



ated within some approximation scheme, e.g.. 
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an additive measure of the deviation of the density matrix of a several-qubit system from 



the "ideal" density matrix of the system 



33, 



3J]. The latter would describe the "ideal" 



dynamics, i.e., the system completely isolated from the environment. For a spin-boson 
model of a single two-state system (a qubit), this measure is calculated explicitly for the 



environment modeled as a bath of harmonic modes j35|, e.g., phonons or photons. We also 
establish that for a several-qubit quantum system, the introduced measure of decoherence 
is approximately additive for the time scales of interest, i.e., for short "gate function" times. 
The fact that this measure can be evaluated by summing up the deviation measures of the 
constituent qubits, which can be, in general, entangled, allows to avoid lengthy, tedious, and 
in most cases, intractable, many-body calculations. This new short-time additivity property 
is reminiscent of the approximate additivity expected for relaxation rates of exponential 
approach to equilibrium at large times, though the two properties are not related. 

Let us briefly outline the commonly accepted scheme for implementing quantum algo- 
rithms in physical systems. The input data are encoded into the quantum states of several 
separated two- level systems (qubits). Then, their evolution is controlled b y a Hamiltonian 
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. Parameters 



consisting of single-qubit operators and of two-qubit interaction terms 
of the Hamiltonian can be varied (controlled) externally to im plem ent a g iven al gori thm. In 



most quantum computer proposals [3£ , 



most quan tum cc 
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3S, 



naiiy to im plem ent a g iven algorii 

3 H Q, H y, HQ, Q, H 
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this control is achieved by changing local electromagnetic fields around the 
qubits. Of course, this ideal model does not include the influence of the environment on the 
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computation, which necessitates quantum error correction. The latter involves inevitably the 
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implementation of non-unitary measurement-type operations 
and cannot be described as Hamiltonian-governed dynamics of a closed system. 

To study the effect of the environment, one needs to choose a suitable model for the 
environmental and its coupling to the system of interest. The accepted approach to evaluate 
environmentally induced decoherence involves a model in which each qubit is coupled to a 
bath of environmental modes SSl- The reduced density matrix of the system, with the bath 
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modes traced out, then describes the time-dependence of the system evolution ja 17|, |8|, |9|, 
0, H 12, Q 14, H 12, 3 3 H 0, Due to the interaction with 

the environment, after each computational cycle the state of the qubits will deviate slightly 
from the ideal state. The deviation accumulates at each cycle, so that large-scale quantum 

e without performing fault-tolerant error correction schemes 



confutation is not possi 
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l62| . These schemes require the environmentally induced 
decoherence of the quantum state in one cycle to be below some threshold. The value of 
threshold, defined for uncorrelated single qubit error rates, was estimated Q to be 
between 10~® and 10~^. 

To study decoherence for a given system, one should first obtain the evolution of its 



density matrix. This can 

Q 0, 



3e done by using various approximations 
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0,0,0,0 

Hamiltonian of the studied system is usually considered to be constant. The most familiar 
are the Markovian-type approximations 00,0,0, used to evaluate approach to the 
thermal state at large times. It has been pointed out recently that these approximations are 
not suitable for quantum computing purposes because they are usually not valid at 
low temperatures_an_d for short cycle times of quantum computation. Several non-Markovian 



During each computational step the 



approaches 
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30| have been developed 



to evaluate the short-time dynamics of open quantum systems. 

When one tries to study decoherence of several-qubit systems, additional difficulties 
should be taken into account. Namely, one has to consider the degree to which noisy 
environments of different qubits are correlated 6^. For example, if all constituent qubits 
are effectively immersed in the same bath, then there is a way to reduce decoherence for 
this group of qubits without error correction algorithms. The reduction of error rate can 
be achieved by encoding the state of one logical qubit in a decoherence free subspace of the 
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states of several physical qubits j66|, |67|, |68|, |69| . Therefore, in a large scale quantum infor- 
mation processor consisting of many thousands of qubits, it is more appropriate to consider 
qubits immersed in distinct baths, because these errors represent the "worst case scenario" 
that necessitates error-correction. 

After obtaining the density matrix p{t) for a single- or few-qubit system evaluated in 
some approximation, we have to compare it to the ideal density matrix p*^*^ (t) corresponding 
to quantum algorithm without environmental influences. It is convenient to define some 
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measure of decoherence to compare with the fault-tolerance criteria 

0, • It is desirable to have the measure nonnegative, and vanishing if and only if the 
system evolves in complete isolation. Since explicit calculations beyond one or very few 
qubits are exceedingly difficult, it would be also useful to find a measure which is additive 
(or at least sub- additive), i.e., the measure of decoherence, D of a composite system will be 
the sum (or not greater than the sum) of the measures of decoherence, Dj, of its subsystems, 

D<J2Dr (1) 

j 

In this work, we identify a measure of decoherence which has these desirable properties. 
In Section InJ we give an overview of different methods for quantifying decoherence. We 
note that in some cases these numerical measures have oscillations on the time scales of 
the internal system frequencies, which do not reflect the nature of decoherence. Therefore, 
in the next Section IIIII we deflne the norm that is subadditive for non-interacting initially 
unentangled qubits and in most cases is a monotonic function of time. To establish stronger 
properties of subadditivity of this norm even for initially entangled qubits we will use a more 
sophisticated diamond norm. It is deflned in Section HVl Relation between these norms and 
conditions for subadditivity for initially entangled qubits are also discussed. In Section FV] 
we consider a speciflc model of a qubit interacting with a bosonic bath of environmental 
modes. For two types of interaction we explicitly obtain norms D and K for one qubit and 
prove asymptotic additivity for short times, t, 



(2) 



This result should apply as a good approximation up to intermediate, inverse-system-energy- 
gap times , and it is consistent with the recent flnding [t^ that at short times decoherence 
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of a trapped-ion quantum computer scales approximately linearly with the number of qubits. 



II. DIFFERENT APPROACHES TO QUANTIFYING DECOHERENCE 



Typically, the total Hamiltonian of an open quantum system interacting with environment 
has the form 

H = Hs + Hb + Hj, (3) 

where Hs is the internal system Hamiltonian, Hb is the Hamiltonian of the environment 
(bath), Hj is the system-bath interaction Hamiltonian. Over gate-function cycles, and 
between them, the terms in H will be considered constant 0, 0] . Therefore, the overall 
density matrix R{t) of the system and bath evolves according to 



R{t) 



-iHt 



R{0)e 



\Ht 



Here and in the following we use the convention h = 1. 
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58, 



(4) 



to be a direct product of the 



Usually the initial density matrix is assumed |56 , 
initial density matrix of the system, p(0), and the thermalized density matrix of the bath, 

e, 

i?(o) = p(o)®e. (5) 



The reduced density matrix of the system is obtained by tracing out the bath modes. 



(6) 



Usually the environment is assumed to be a large macroscopic system in thermal equilib- 
rium at temperature T. Interaction with it leads to thermalization of the quantum system, 
so that the reduced density matrix at large times approaches 



p 



Tr 



s le" 



as t 



oo. 



(7) 



where /3 = 1/kBT. Markovian type approximations which are used to quantify this process 
yield exponential decay of the density matrix elements in the energy basis of the Hamiltonian 
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(8) 



p^^{t) Qc e-*/^"™ (n 7^ m). 



(9) 
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The shortest times among T„„ and T„^m are identified as characteristic times of thermahza- 
tion Ti and decoherence T2, respectively. Thermahzation of the quantum system requires the 
energy exchange with the environment while decoherence can include other faster processes 
without energy exchange. Therefore, it is commonly believed js^ that T2 is much shorter 
than Ti for low-temparature, well isolated from the environment systems appropriate for 
quantum computing realizations. The shorter time T2 can be compared with Tg needed for 



elementary quantum gate functions 71[. The ratio Tg/T2 must be small in order to satisfy 
fault-tolerance criteria [gsI . 6^ . However, since T2 is a large-time asymptotic property, other 
measures, representative of the short-time, t <tiT2, decoherence properties, are preferred for 
actual numerical evaluations jisl . 

The exponential behavior of the density matrix elements inherent to the Markovian ap- 
proximation is valid on the time scale which is much larger than the internal inverse-system- 
energy-gap times of the quantum system. To measure effects of decoherence on these time 
scales one can use the entropy jl^, 

5(t) = -Tr(plnp), (10) 

or the idempotency defect, called also the first order entropy 

sit) = 1 - Tr (/) . (11) 

Both expressions are zero only if the quantum system density operator is a projector p(0) = 
\ip){(p\. Any deviation from a pure state leads to the increase in the value of both measures. 
Expressions (fTTH ITT|l provide numerical measure of the system "purity" which does not rely 
on preferred basis. However, entropy measures do not distinguish different pure states. 

The next step in analyzing the effect of the interaction with the environment is to define 
the "ideal" (without interaction) density operator evolution according to 

p»(t) = e-'^^V(0)e'''^*. (12) 

One of the measures that characterizes decoherence in term of the difference between the 
"real" evolution, p(t), and "ideal" one, p^^\ is the fidelity 

F{t)=TTs[p^Ht)p{t)]. (13) 

It is particularly useful when p'^'^\t) remains a projection operator (pure state) for all times 
t > 0, because it then attains the maximum value of 1 only when pit) = p^^\t). 



An alternative way to quantify the effect of the interaction with the bath j32], is to 
consider the deviation, cr{t), 

a(t)=p(t)-p«(t), (14) 

of the reduced density matrix from the ideal one. As numerical measures of decoherence one 
can use the operator norni ||cr|K or trace norm ||o"||xr- Both norms are standard in in the 
theory of linear operators [77|. For an arbitrary linear operator, A, these norms are defined 
as follows, 

piL-upfM^y^ (15) 

\\A\\^ = Tt Va^. (16) 

For a finite-dimensional Hermitian deviation operator ()14j) these definitions are equivalent 

to 

\\o-\\x = inax |Aj| , (17) 



l^llTr 



where Aj are the eigenvalues of a. 

In the simplest case of a two- level system (qubit), the two norms are proportional to each 
other and given by 



^IIa = Vl^lll^ + 1^121^ = \ ll^llxr- (19) 



The deviation norm WcW^ has the minimal value, 0, only for p(t) = p^'^\t) without any 
additional conditions. 

Note that the measures ||cr||^ and \\cr\\^ are not only functions of time, but also depend 
on the initial density operator p(0). Due to decoherence, they will deviate from zero for 
t > 0. However, their time- dependence will also contain oscillations at the frequencies of the 
internal system dynamics, as will be illustrated below for ||cr||;^; see 32[|34| and Fig. 1. In 
the next section we define the maximal deviation norm, D{t), which is typically a monotonic 
function of time. 
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III. THE MAXIMAL DEVIATION NORM 



To characterize decoherence for an arbitrary initial state, pure or mixed, we propose to 
use the maximal norm, which is determined as an operator norm maximized over all the 
possible initial density matrices. For instance, we can define 

D(t)=supf l|a(t,p(0))l|,y (20) 
p(o) V / 

One can show that < D{t) < 1. This measure of decoherence will typically increase 
monotonically from zero at t = 0, saturating at large times at a value D{oo) < 1. The 
definition of the maximal decoherence measure D{t) looks rather complicated for a general 
multiqubit system. However, we will show that it can be evaluated in closed form for short 
times, appropriate for quantum computing, for a single-qubit (two-state) system. We then 
establish an approximate additivity that allows us to estimate D{t) for several-qubit systems 
as well. 

In the superoperator notation the evolution of the reduced density operator of the syst em 
dni) and the one for the ideal density matrix (|T^ can be formally expressed 0, 0, 0| in 
the following way 

p(t)=T(t)p(0), (21) 

p«(t)=r«(t)p(0), (22) 
where T, T*-*) are linear superoperators. In this notation the deviation can be expressed as 

a(t)= [T(t)-T«(t)]p(0). (23) 

The initial density matrix can always be written in the following form, 

P(0) = 5^P,|^,)(^,|, (24) 
j 

where ^jPj = 1 and < < 1. Here the set of the wavefunctions is not assumed to 
have any orthogonality properties. Then, we get 

^ {t.pm = J2p^ [m - r«(t)] iv^,) . (25) 
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The deviation norm can thus be bounded, 

lk(t,p(O))||A<||[T(t)-T«(t)]|0)(0|||,. (26) 

Here |0) is defined according to 

II [T - T«] 10) (01 II, = max || [T - T«] |^,)(^,|||, • 

It transpires that for any initial density operator which is a statistical mixture, one can 
always find a density operator which is pure-state, |0)(0|, such that ||cr(t, p(0))||a < 
||cr(t, |0)(0|)||a. Therefore, evaluation of the supremum over the initial density operators 
in order to find D{t), see (j2DI), one can done over only pure-state density operators. 

Let us consider strategies of evaluation of D{t) for a single qubit. We can parameterize 
p(0) as 

p(0) = f/(^ " ){/', (27) 

where < P < 1, and U is an arbitrary 2x2 unitary matrix, 

/ e^("+T)cos^ e^("-^)sin^ \ 
I _e*(7-") sin^ e-*("+^) cos^ / 

Then, one should find a supremum of the norm of deviation p7p over all the possible real 
parameters P, a, 7 and 6. As shown above, it suffices to consider the density operator in 
the form of a projector and put P = 1. Thus, one should search for the maximum over the 
remaining three real parameters a, 7 and 6. 

Another parametrization of the pure-state density operators, p(0) = \4>){4>\, is to express 
an arbitrary wave function \(f)) = + some convenient ortonormal basis \j), 

where j = 1, . . . , N. For a two-level system, 

where the four real parameters aj, bj, with j = 1, 2 satisfy af + bf + a2 + b"^ = 1, so that the 
maximization is again over three independent real numbers. In the following sections, we 
will consider examples of evaluation of D{t) for single-qubit systems. 

In quantum computing, the error rates can b e sig nificantly reduced by using several phys- 



ical qubits to encode each logical qubit 



67 



68, 



69|. Therefore, even before active quantum 



error correction is incorporated 20|, |2l|, |22, |23|, |2J, |25|, |26|, |27|, evaluation of decoherence of 
several qubits is an important, but formidable task. Consider a system consisting of two 
initially unentangled subsystems Si and 5*2, with decoherence norms Ds^ and Ds2i respec- 
tively. We denote the density matrix of the full system as psiS2 ^"^^ deviation as cr^^^^i 
and use a similar notation with subscripts 5*1 and 5*2 for the two subsystems. If the evolution 
of system is governed by the "noninteracting" Hamiltonian of the form Hs^S2 = -^^i + Hs2-i 
where the terms Hs^, act only on variables of the system 5*1, 5*2, respectively, then the 
overall norm D 3-^^32 can be bounded by the sum of the norms Ds^ and Ds2, 



Ds,S2 = sup IksisJU = sup \\ps,s2 - pSsJIa 

p(0) p(0) 



(30) 



sup \\psfPs2 - Psl'^Psl II A = sup \\asfPs2 + Psl^(^S2 II a 



P(0) 



P(0) 



Since the operator norm obeys the triangle inequality jl^l 

1 1 ^5i®P52 + Psl J I A < 1 1 ^5i®P5j I A + 1 1 P£ ®^52 1 1 A , 

we can estimate the last expression as 



(31) 



S^S2 = sup lk5i^P52+ P5>^52 II A < SUp H^T^^^P^JU + SUp UpljJ^a^JI a- 
p(0) p(0) p(0) 



(32) 



Each eigenvalue of the tensor product of two linear operators is formed as a pairwise product 
of eigenvalues of the two operators. Therefore, the operator norm of the tensor product of 
two operators is equal to product of their operator norms. 



\Ai^A2\ 



I^IIaP2||a. 



(33) 



We use this property and the fact that the eigenvalues of density matrices ps^ , ps2 are in 
[0, 1] to derive the estimate 



Ds,s2< sup \\asjx+ sup \\asjx = D^^ + D^^. 

PSi(O) PS2(0) 



(34) 



In general, initially unentangled qubits will remain nearly unentangled for short times, 
because they didn't have enough time to interact. Therefore, the inequality 



(35) 
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is expected to provide a good approximate estimate for the norm of a multiqubit system D 
in terms of the norms Dq calculated in the space of each individual qubit, i.e., the measures 
of decoherence for the individual qubits can be considered approximately additive. For large 
times, the separate measures become of order 1, so such a bound is not useful. Instead, the 
rates of approach of various quantities to their asymptotic values are approximately additive 
in some cases. 

In the rest of this work, we focus on the short-time and adiabatic (i.e., no energy ex- 
change with the bath regimes, and establish a much stronger property: We prove the 
approximate additivity for the initially entangled qubits whose dynamics is governed by 

H = Y,Hq = Y,{Hsq + Hb, + Hiq) , (36) 

where Hsq is the Hamiltonian of the gth qubit itself, is the Hamiltonian of the envi- 
ronment of the gth qubit, and Hiq is corresponding qubit-environment interaction. For this 



purpose 
norm l60 



in the next section we consider a more complicated (for actual evaluation) diamond 



61 



62|, K{t), as an auxiliary quantity used to establish the additivity of the more 



easily calculable operator norm D{t). 



IV. THE DIAMOND NORM 



The establishment of the upper-bound estimate for the maximal deviation norm of a mul- 
tiqubit system, involves several derivations. We bound this norm by the recently introduced 
(in the contexts of quantum computing) 6^ 61, 131 diamond norm, K{t). Actually, for 



single qubits, in several models the diamond norm can be expressed via the corresponding 
maximal deviation norm. At the same time, the diamond norm for the whole quantum 
system is bounded by sum of the norms of the constituent qubits by using a specific sta- 
)ilit y p roperty of the diamond norm. The use of diamond norm was proposed by Kitaev 



6G 



yjp 



6^ 



Kit) = \\T - T»|U = sup ||{[r - T«]«/}f?||Tr. (37) 



e 



The superoperators T, T*^*) characterize the actual and ideal evolutions according to (j2Ill22I), 
I is the identity superoperator on a Hilbert space G whose dimension is the same as that 
of the corresponding space of the superoperators T and T*^*\ and q is an arbitrary density 
operator in the product space of twice the number of qubits. 
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The diamond norm has an important stabihty property, proved in 

|60|,|6lll6^, 



in leO 



||Ti«T2||o= IIT1IUIIT2II0, (38) 

which should be compared with ()33p. Note that ()38p is a property of the evolution super- 
operators rather than that of the density operators: The importance of this difference will 
become obvious shortly. 

Consider again a composite system consisting of the two subsystems Si, S2, with the 
noninteracting Hamiltonian 

HsiS2 = + Hs2- (39) 
The evolution superoperator of the system will be 

Ts,s,=Ts,^Ts,, (40) 

and the ideal one 

T^s.=T^<- (41) 
The diamond measure for the system can be expressed as 

By using the stability property we get 



\\T^'!\\o\\Ts, -Ti^\\o= \\Ts, - T«|U + \\Ts, - = K^^ + K^.^. (43) 

The approximate inequality 



K<Y,K,, (44) 



for the diamond norm K{t) has thus the same form as for the norm -D(t), (jH^ . Let us em- 
phasize that both relations apply assuming that for short times the subsystem interactions, 
directly with each other or via their coupling to the bath modes, have had no significant 
effect. However, there is an important difference in that the relation for D{t) further requires 
the subsystems to be initially unentangled. This restriction does not apply for the relation 
derived for K{t). This property is particularly useful for quantum computing, the power of 
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which is based on qubit entanglement. However, even in the simplest case of the diamond 
norm of one qubit, the calculations are extremely cumbersome. Therefore, the measure D(t) 
is preferrable for actual calculations. 

The two deviation-operator norms considered Section O are related by the following 
inequality 

M,<^M^,<1. (45) 
Here the left-hand side follows from 

Tra = J]Aj = 0. (46) 
j 

It follows that the ith eigenvalue of the deviation operator a that has the maximum absolute 
value, Xi = Amax, can be expressed as 



Xe = -J2Xr (47) 



Therefore, we have 



\WL = \ (2|A.i) <\[\M + Y. i^^i) = \ i^^i j = 2 ii^iItv • (48) 

The right-hand side of ()45|) then also follows, because any density matrix has trace norm 1, 

Iklk = \\P - P^'^hr < \\phr + ll/^lliY = 2. (49) 
From the relation PU)) it follows that 

K{t) < 2. (50) 
By taking supremum of the both sides of the relation (|lHj) we get 

D{t) = sup < i sup llall^^ < ^K{t), (51) 

p(0) ^ p(0) ^ 

where the last step involves technical derivation details not reproduced here. In fact, in 
the following sections we show that for a single qubit, calculations within selected models 
actually give 

D{t) = ^K{t). (52) 
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Since D is generally bounded by (or equal to) it follows that the multiqubit norm D 

is approximately bounded from above by the sum of the single-qubit norms even for the 
initially entangled qubits, 

D{t)<\K{t)<\Y.^<, = Y,D,, (53) 

where q labels the qubits. 



V. DECOHERENCE IN THE SHORT-TIME APPROXIMATION 

Typically the environment, a large macroscopic system, is modelled by a bath of an 
infinite number of modes. Each mode is represented by its own Hamiltonian M^, 



Hb = Y^ Mk. (54) 



k 

The interaction with the bath is often described by the coupling of its modes to Hermitian 
operator of the quantum system. 



Hi = AsJ2'Jk. (55) 



k 

For a bosonic-mode heat bath |78|1 we take 

Mfc = uJkala^, Jk = g^al + gla^. (56) 

Here Uk are the bath mode frequencies, a^, a[ are the bosonic annihilation and creation 
operators, and gk are the coupling constants. Two eigenbases of the operators Hs and As 
are 

Hs\n) = En\n), Ash) = X^h). (57) 

At the the initial time t = the total density matrix of the system and bath is a direct 
product -R(O) = p(0)®B of the initial density matrix of the system p(0) and the density 
matrix of the bath G. The latter is a product O = 6*1® ^2 ■ ■ ■ of the bath modes density 
matrices Ok- Each bath mode k is assumed to be thermalized, 

g-/3Mfc 
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In the short-time approximation the exponentials in ^ representing the time evolu- 
tion of the total density matrix R{t) are approximated as 

The matrix elements of the reduced density matrix p{t) in the free energy basis can be 
expressed as 

After cumbersome calculations utilizing the fact that the trace over the bath modes 
can be carried out separately for each mode, for the bosonic bath case one obtains 



B2(t)(A^„A.)2/4+iC{t)(A2-A?) 



Pmn \ 

(61) 



where 

2 

\gk\ 2 PuJk 



5^(t) = 8$:^sin^^^coth^, (62) 



k -"fc 

|2 



C{t) = y2 —^{ukt - sinwfct). (63) 

k 

Here the Roman-labeled states, are the eigenstates of Hs corresponding to the eigenvalues 
Ei, with i = m,n,p,q. The Greek-labeled states, \(), are the eigenstates of with the 
eigenvalues A^, where ( = fijU. Details, more general expressions, and additional discussion 
can be found in 



VI. THE SPIN-BOSON MODEL 

Let us consider a system which is a spin-1/2 particle in an applied magnetic field, inter- 
acting with the boson bath. In this case the system Hamiltonian is 

Hs = (64) 

and interaction can be chosen as As = (Tx, where ax and are the Pauli matrices, and 
i7 > is the energy gap between the ground (up, |1) =|t)) and excited (down, |2) =|i)) 
states of the qubit. The eigenstates of ax will be denoted by 

ax\±) = ±\±), (65) 

15 



where 



l±)^^(|l>±|2» 



The dynamics of the system can be obtained in closed form as 



,(t)= J2 {m\fi){fi\p){qW){u\n)p,,{0)e^(-- 



m)5mp5nq^t-B'^{t)(l-S^^) 



(66) 



(67) 



/j,,!^ — ih 



where 6mn is Kronecker symboL 

Note that this resuh depends only on the spectral function B'^{t), defined in because 



(68) 



This function is obtained by integration over the bath mode frequencies. When the sum- 
mation in (|62|) is converted to integration in the limit of infinite number of the bath modes 
BSQ, we get 



B^{t) = 8 J diuN{iu)\g{iu)\''iu-^sm^^cot\i^, 



(69) 



where N{uj) is the density of states. In many realistic models of the bath, the density of 
states increases as a power of u for small frequencies and has a cutoff, Uc, at large frequencies 
(Debye cutoff in the case of a phonon bath). Therefore, approximately setting 

N{uj)\g{uj)\'^ (xuj''exp{-uj/ujc) (70) 

can yield a good qualitative estimate of the relaxation behavior P, . For a popular case 
of Ohmic dissipation n = 1 and the function B'^{t) has the initial stage of quadratic 
growth, intermediate region of logarithmic growth, and linear-in-t large-time behavior. 
Evaluation of (|H7jl yields the following expressions. 



P22{t) 



1 + e 



-BHt) 



P22(0) 



1 - e^^'W 



Pii(O) 



(71) 



P2l(t) 



-iQt 



1 + e 



P2l(0) 



1 - e 



Pl2(0) 



Deviation operator a (t) is defined by 



0-22 (i) 



1-e 



[Pll(0)-P22(0)] 



(72) 



(73) 
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a2l(t) 



-BHt) 



[pi2(0)-e-^V2i(0)] 



(74) 



With pi2(0) = |pi2(0)|e*'^, we get 



I [1 - e-^^«] {[pn(0) - P22(0)]^ + 4 ^0)]' sin2[(fi/2)t + (75) 



In Fig. 1, we show schematically the behavior of ||cr(t)||^ for three representative choices of 
the initial density matrix p(0). Generally, the norm ||cr(t)||^ increases with time, reflecting 
the decoherence of the system. However, oscillations at the system's internal frequency fl 
are superimposed, as seen explicitly in (f7H|) . Thus, the decohering effect of the bath is better 
quantified by the maximal operator norm, D{t). Explicit calculations yield the result, shown 
in Fig. 1, 



1 - e-^'W 



(76) 



which is indeed a monotonically increasing function of time. 

Let us now consider the diamond norm K{t) for one qubit. We will later use the result 
to estimate the norm D{t) for the multi-qubit case. To find K{t) for a two-level system one 
has to deal with the 4x4 density matrix gjk,im, where j, k,l,m = 1, 2. To evaluate 



r^{(r-r«)«/}f?. 



(77) 



one assumes that T — T^*) acts in the subspace labeled by the indices j, I, see fl73l74|) . while 
the subspace labeled by the remaining pair of indexes, k, m, is unaffected. The resulting 
expression are 



^2fc,2m 



[f?lfe,lm ~ Q2k. 



2m 



(78) 



^2k,lm — 2 



\_Qlk,2n 



-int 



Q2k,lm\ 



lfc,lm 



2k,2m, 



E- 



lfc,2m 



y 



2m,lfc5 



k,m = 1, 2. 



The maximal trace norm of the 4x4 matrix E is calculated by considering the pure-state 
density matrices g = similarly to the consideration in Section ITTTl see (|23II29|) . Here 

10) can be expressed in terms of the basis states. 



(79) 



m,n=l,2 
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t 



FIG. 1: The maximal deviation norm D vs. time for the spin-boson model with Ohmic bath. The 
non- monotonia curves illustrate the behaviour of the norm for several initial choices of 

density operator /9(0). 

The constants a^^, are normalized real amplitudes, such that Ylimn=i 2('^nm + ^nm) = 1- 
The eigenvalues of 17, see (I78p . are = and 



<^3,4 = ± 



"1 _e-^'W" 








2 





(011021 + 012022 + ^11^21 + &I2&22) COS ( ^ ] 



- (021^11 + 022&12 - 011621 - 012622) sin ( 



2^ 1/2 



(80) 



The maximal value corresponds to the square root in the expression (jHn|) equal to 1, and 
thus the diamond norm is 

ir(t) = 1-e-^'W. (81) 
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The above calculation establishes that K{t) = 2D{t) for the spin-boson model at short 
times, and gives the upper bound on the multiqubit norm D, 

^^E^^ = ^E[i--""^^1- (82) 

For short times, one can also establish a lower bound on D(t). Consider a specific initial 
state with all the Q qubits excited, p(0) = (|2) (2|)i® . . . '»(|2) (2|)q. Then according to 
dZHZa), pit) = pi(t)0 . . . 0pQ{t), where 

1 



Pgit) 



1 + e-^^ 



)(2|}. 



1 - e-"« 

The right-bottom matrix element of the (diagonal) deviation operator 
can be expanded, for small times, as 



(83) 



g 



(84) 



q \ 1 / 



(85) 



because B{t) vanishes for t — 0. The largest eigenvalue of a{t) cannot be smaller than 
c"2<3,2<3(^)- It follows that 



g \ g / g \ g / 

where we used flTUj) for short times. 

By combining the upper and lower bounds, we get the final result for short times, 

g \ g J 



(86) 



(87) 



VII. SPIN MODEL FOR PURE DECOHERENCE 



Let us consider a two-level system interacting with the bosonic bath of the environmental 
modes adiabatically, i.e., without exchange of energy with the bath modes 79|. While 
energy exchange processes are needed for thermalization, and also contribute to decoherence, 
additional, "pure" decoherence processes are possible and are expected to be important at 
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low temperatures and short to intermediate times, appropriate for quantum computing 
designs [l8|. The spin-model Hamiltonians Hs and Hb will still be assumed of the form 

Hs = -%z, (88) 



Hb = '^ujkalak. (89) 

k 

However, the interaction term will be now of the form commuting with the system's energy, 
thus making the latter conserved, 

Hi = cT.J2 {9k4 + ■ (90) 

k 

Since Hs commutes with Hj, instead of the approximate formula (j^^ we have the exact 
factorization. 

Therefore, the analytical expression for the reduced density operator (jfilll is exact and valid 
for all times . For the two-level case 

Pii(t)=Pii(0), (92) 



P22(t)=P22(0), (93) 

Pi2(t) = Pi2(0)e^™-^'«, (94) 

where the spectral function B{t) was defined in (jU^ . The diagonal density matrix elements 
remain constant because there is no energy exchange of the system with the bath. However, 
there is pure (adiabatic) decoherence manifested by the decay of the off-diagonal elements, 
characterized by the spectral function B[t) The corresponding deviation matrix is 

expressed as follows, 

fTll(t) =a22(t) =0, (95) 

ai2 = Pi2(0)e*™-^'W. (96) 
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The operator norm of cr is 



1 - e-^'W 



IP12I 



(97) 



Since both diagonal elements of the density matrix and its eigenvalues are in [0, 1] it follows 
that the absolute value of the off-diagonal element of any two-level-system density matrix 
less than 1/2, 

(98) 



IP12I < \. 



Thus, 



D(t)=sup(||a(t,p(0))|U) 



p(o) 



1 - e-^'W 



One can also evaluate 



r = {(r-T«)«/}f?, 



(99) 



(100) 



jk,jm — 0, 



;ioi) 



For g 



Since 



nk,2m — ^2mAk 



^ (()\jnt-B'^{t) 

Qlk,2m[^)e 



, with 10) as in (ffn|) . the eigenvalues of E are <;i = ^"2 = and 



%4 = ± 



1 - e-^'W 



,^?ll,ll + ^'12,12) (^?21,21 + ^•22,22)- 



one can show that 



with 



-^3,4 — i 



Tr^ = 1, 

1 - e-^^W" 



x(l — x) 



X — ^?ll,ll + ^"12, 12 

satisfying < x < 1. The diamond norm (jHTjl thus follows, 



K{t) = sup <{ 2 

0<x<l 



1 - e' 



a/x(1 — x) [ = 1 — e 



(102) 



(103) 



(104) 



(105) 



(106) 



(107) 



It is instructive to compare ()99I107|) and ()76I8H) . The results are identical despite the fact 
that the interaction terms and density matrix time- dependence ()71l72l92l94j) are different. 



21 



As in the case of the short-time approximation, we get the upper bound for the multiqubit 
norm D{t), 

To estabhsh the lower bound on D{t), we consider a specific initial state with p(0) = 
|\E')(\1'|, which is a superposition of the state corresponding to all the qubits in their ground 
states and that of all qubits in their excited states, 

1 



Then according to (j92l94jl . 



|vI/) = — (|1...1) + |2...2)). 



Pl,l{t) = p2Q,2Q{t) = ^, 



(108) 



(109) 



(110) 



Pi,2Q (^) = ^ exp -i^Qqt-^B^g [t) 

q q 

The only non-zero matrix elements of the deviation operator are the right-top and left- 
bottom matrix elements. 



(Ti_2Q(t) 



1 _e-^.^9^ 



;iii) 



For short times, the absolute value of c'"i 2Q(t) can be expressed as 



q \ q / 



(112) 



where the first term gives the largest eigenvalue of cr{t). It follows that 



q \ q / q \ q / 



(113) 



where we used (I99|) for short times. Finally, we get the same result (j87|) for the approximate 
additivity of D for short times, for the present model of adiabatic decoherence. 



q \ q / 



(114) 



In summary, we introduced the maximal operator norm suitable for evaluation of deco- 
herence for quantum system immersed in a noisy environment. The new maximal operator 
norm was evaluated for spin models with two types of bosonic bath interaction. We es- 
tablished both general and model specific subadditivity and additivity properties of this 
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measure of decoherence for multi-qubit system at short times. The latter property allows 
evaluation of decoherence for complex systems in the regime of interest for quantum com- 
puting applications. 

This research was supported by the National Security Agency and Advanced Research 
and Development Activity under Army Research Office contract DAAD-19-02-1-0035, and 
by the National Science Foundation, grant DMR-0121146. 
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